NASA Technical Memorandum 85912 

NASA-TM-85912 19860015695 


Combining and Connecting Linear, 
Multi-Input, Multi-Output Subsystem 
Models 


Eugene L. Duke 


LIBRARY COPY 

Sr'P 2 3 1986 

M1QLEY RESEARCH CENTER 
LIBRARY, NASA 
HAMPTON, VIRGINIA 


IW\SA 

National Aeronautics and 
Space Administration 





NASA Technical Memorandum 85912 


Combining and Connecting Linear, 
Multi-Input, Multi-Output Subsystem 
Models 


Eugene L Duke 

Ames Research Center, Dryden Flight Research Facility, Edwards, C&lifornia 


NASA 

National Aeronautics and 
Space Administration 
Ames Research Center 
Dryden Flight Research Facility 
Edwards, California 93523 


ft? 






SUMMARY 


This report provides the mathematical background for combining and connecting 
linear, multi-input, multi-output subsystem models into an overall systems model. 
Several examples of subsystem configurations are examined in detail. The report 
also contains a description of a MATRIX x command file to aid in the process of com- 
bining and connecting these subsystem models. 


INTRODUCTION 


To analyze and design models for complex systems, it is essential to combine 
subsystem models into an overall system model. This report provides the mathemati- 
cal background for combining and connecting linear, multi-input, multi-output sub- 
system models. Throughout this report, systems are assumed to have a standard form. 

x = Ax + Bu 
£ = Cx + Du 

Using this system model, the mathematics that combine and connect two systems at 
a time are described. Six representative combinations are presented: parallel 

systems (fig. 1), summation of systems (fig. 2), concatenation of systems (fig. 3), 
systems with a common input (fig. 4), summation of systems with a common input 
(fig. 5), and feedback systems (fig. 6). The matrices representing the combined 
system are presented in terms of the matrices for the subsystems. 

The mathematics used in the development of the sample systems are presented in 
appendix A in a MATRIX x (ref. 1) command file format. A listing of the command file 
and instructions for the use of the file are presented in this report. MATRIX x is a 

general purpose, interactive program for matrix manipulation, control law analysis 
and design, and parameter identification.* 


NOMENCLATURE 


n 

dimension 

of 

state vector, x 

m 

dimension 

of 

control vector, u 

r 

dimension 

of 

observation vector, % 

s 

Laplace variable 


*MATRIX x is a proprietary product of Integrated Systems Inc., Palo Alto, 
California. 


Matrices : 


A state matrix 

B control matrix 

C observation matrix 

D feedforward matrix 

F input selection matrix 

H output selection matrix 

I identity matrix 

N (I + D-j H 2 D 2 F 2 ) -1 

0 matrix of zeros 
Vectors : 

u control vector 

w selected input vector 

x state vector 

observation vector 
z selected output vector 

Subscripts : 

1 system one 

2 system two 

upper left submatrix of a partitioned matrix 
2_ upper right submatrix of a partitioned matrix 
3^ lower left submatrix of a partitioned matrix 
4 lower right submatrix of a partitioned matrix 
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MATHEMATICAL PRELIMINARIES 


Matrices, which are a convenient method of formulating simultaneous linear 
equations, are extremely useful in representing the linear models of dynamic 
systems. Just as a system of linear scalar equations can be transformed into a 
matrix equation, a system of linear matrix equations can be transformed into a 
matrix equation using partitioned matrices. Thus the process for converting the 
scalar equations 


ail X, + . . . + a 1n x n = b 1 


a n i x 1 + * • • + a nn x n b n 


into the matrix equation 


a 11 • • • a 1 n - 


“ X 1 


1 

• tr 

— k 

1 

• 


• 

= 

• 

_ a n-| • • • a n n _ 


_ x n_ 


_*>n_ 


is identical to that used for converting the matrix equations 


A 1-l X 1 + 


+ A 1 n £n = kl 


A n i X 1 + 


+ Ay 


2Sn bjj 


into the partitioned matrix equation 


A 1 1 1 • • • 1 A 1 n 


X 1 


b i 

| 1 

1 • 1 





I * t 

1 * 1 

_ A nT | * • * | A n n _ 


• 


• 


_ x n_ 




Partitioned matrices are manipulated in the same way as normal matrices except 
that the matrices and partitioned submatrices must be conformable for the operation 
being performed. Consider the addition of matrices A and B, each partitioned into 
four submatrices. 


A 


~ A J_ | a 2 
a 3 I a 4 
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B 


_B 1_ I B 2~ 

_ B 3 j" b 4_ 

If A and B have the same order and the corresponding submatrices (for example, A-| 
and B-|) are of the same order, the sum may be expressed in terms of the partitioned 
submatrices. 


_A 1_ + B J_ | a 2 + b 2~ 

A + B = B - H - A ~ “ £•£ 

The same principles hold for multiplication. For example, if matrices A and B 
are conformable for multiplication and all appropriate partitioned submatrices are 
conformable for multiplication, then the product can be represented as 


AB 


A 1 B 1 + A ^ B 3_ I a 2 b j 2 + a 2_ b 4 
a 3 b 4 + a 4 b 2 "j* A 2 B 2 + A ! B 4 


DEVELOPMENT OF COMBINED SYSTEM MODELS 


In this section, the models for six representative subsystem combinations are 
developed. These subsystem combinations are for two independent (parallel) systems 
(fig. 1), the summation of two systems (fig. 2), the concatenation of two systems 
(fig. 3), two systems with a common input (fig. 4), the summation of two systems 
with a common input (fig. 5), and a feedback system (fig. 6). Each subsystem con- 
sidered is assumed to be represented by linear system equations of the form 

x = Ax + Bu 

y = Cx + Du 


Hi 



=> * 


x 2 = A 2 x 2 + B 2 u 2 
y 2 =C2X 2 +D 2 u 2 


& 


Figure 1. Parallel systems. 
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Figure 5 . Summation of systems with a common input. 


5 

















Figure 6. Feedback system. 

The state vector x is assumed to be an n x 1 vector. The control vector u and the 

observation vector y are assumed to be m * 1 and r * 1, respectively. 

For the systems in which subsystems are connected, the concept of a selection 
matrix is included. For notational convenience, the concept of an input and 

an output selection matrix is used. These matrices are denoted F and H, respec- 

tively, and are related to the subsystem control and observation vectors by the 
following definitions. 


u = Fw 


= Hy 


The vector w represents a generalized input into the system being considered. The 
vector _z serves the same purpose for output. 


For example, an output selection matrix that is not the identity matrix could be 
used as a feedback controller model that is developed independently of the plant 
model. If the feedback control vector, u, of the controller model were 



and the observation vector £ of the plant model were 



a 


the selection matrix would be used to select the last three elements of tjie plant 
observation vector and reorder them to conform to the order of the feedback control 
vector for the controller. For this example, the selection matrix H would be 



"0 

1 

0 

0 " 

H = 

0 

0 

0 

1 


_0 

0 

1 

0 _ 


6 





If scaling of output parameters is desired, the elements of the selection matrix can 
be used to accommodate that scaling. The use of the input selection matrix F is 
identical to that described previously for the output selection matrix. 


Parallel Systems 


Two independent subsystems having no interconnections can be directly combined 
into a single system model simply by defining the new control, state, and obser- 
vation vectors and then, by inspection, constructing the partitioned matrices of the 
resultant combined system. Figure 7 illustrates two parallel systems governed by 
the following equations. 


ill 

£1 


—2 


—2 


+ B-] u-| 

c l2Ll + D liil/ 


System 1 


*2^2 + “A 


C 2—2 + D 2 ± 2 ) 


f System 2 



h 


h 


O 


o 

*1 


o 


*1 


Vi i 0 ' 

i 


*1 

"o” 

-2 


o a 2 


—2 

+ 

0 :b 2 f 2 


—2 


h 


" C 1 i 0 

1 


Si. 

_1_ 

D i F i i 0 ' 

i 


«1 

12 


0 C 2 


XI 

1 

T 

0 ! D 2 F 2 

1 


—2 


Figure 7. State representation of 
combined parallel systems. 
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The control vectors to these systems are defined in terms of the input vectors w-| 
and w , using the input selection matrices F-j and F 2 . 

Hi = F iHl System 1 

H-2 = F 2 — 2 S y stem 2 

The control, state, and observation vectors are defined to include all elements 
contained in the subsystem models. 


u = 



x 



)L 



If the system state equations are expressed in terms of the input vectors and the 
input selection matrices, the equations for the system become 

= A-| x i + B-jF -]W-| System 1 

H 2 = A 2 x 2 + B 2 F 2 w_ 2 System 2 

Using the previous definitions, the new state equation is represented by the 
following equations. 


Hi 


-Ai | 0 - 


HI 


-BiFt , 0 - 

"Hi ” 

• 

= 

j 


— — 

+ 

1 

— 

2 j 


° j a 2 _ 


2_ 


o 1 b 2 f 2 

2_ 

"Z1~ 


~Ct I o - 


"Hi” 


‘DlF-, | 0 - 

"Hi” 

— 

= 

j 


— 

+ 

-j- 

— 

J2_ 


_° I C 2 _ 


2_ 


0 I D2F 2 

_ — 2_ 


The zero submatrices are of an order appropriate for the application. 


Summation of Systems 

Figure 8 illustrates the summation of two systems having independent input vec- 
tors. As shown in the figure, the systems are preceded by an input selection matrix 
and followed by an output selection matrix. The combined system model is developed 
by defining the control, state, and output vectors for the total system; expressing 
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the subsystem controls in terms of the input vectors and input selection matrices ; 
and then expressing the output vector of the total system in terms of the selection 
and observation equation matrices for both subsystems. 



Figure 8. State representation for 
summation of systems . 


The system equations for the two subsystems are given as 




x-| = A-|X1 + B^ui 
Y1 = c-|xi + d-|ui 


System 1 


-2 A 2-2 + B 2-2 ) 
*2 ‘ C 2— 2 + ° 2“2 { 


System 2 


The control vectors for the two systems are defined in forms of input vectors and 
input selection matrices. 


ui = F -|W i System 1 


-2 = F 2-2 System 2 

Substituting these expressions into the state equations for both systems yields 
equations in terms of the system inputs, the system selection matrices, the original 
state vector, and the state and control matrices. 

xi = A-|Xi + B-|F-|wi System 1 
-2 = A 2-2 + B 2 F 2-2 System 2 
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The output of the total system, y, is the sum of the selected outputs of system 
1 and system 2 , z_-\ and z^, respectively. 

y = z.i + z 2 

The selected vectors are defined as 

21 = 

— 4* H-jD -j F -j 

2* = 


System 1 


System 2 


2 ” 2*^2 

‘ H 2 C 2— 2 + h 2» 2 F 2 » 2 J 

Thus the output vector for the total system can be written 
y - HiCi£i + H2C 2 x 2 + H1D1F1W1 + H 2 D2F 2 w 2 

The control vector u and state vector x in the total system are defined as all 
elements that are to be added in the subsystems. This definition yields the 
following partitioned vectors. 


u = 


"1 

- 2 J 

*1 


L* 2 J 


The state equation can be written in terms of these vectors by inspection. 


” • “ 
211 


"at 

1 0 “ 

"211" 


r b i f i 

0 1 

"nr 

• 

_ 

— — 

i 

— 

+ 

— •) 

— — 

L~2_ 


0 

j a 2 

2_ 


L 0 

| b 2 f 2J 

2_ 


The observation equation for the total system can be written using the previous 
definitions for the control and state vectors for the total system and the expanded 
observation equation. 


X = [ H lCl | H 2 C 2 ] 


211 

* 2 ~J 


+ [h^Ft I H 2 D 2 F 2 ] 


, w , 

L-2J 
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Concatenation of Systems 


Two systems concatenated in a series are illustrated in figure 9. The first 
system contains an input selection and output selection matrix. The total system 
representation is determined by defining the state, control, and observation vectors 
for the total system, defining the output selection equation, and defining the 
observation equation from the first system to establish new system equations for the 
second system. 



o 


o 


1 

o 

< 

1 


-1 

4 - 

' B lFl 

0 

—2 


_b 2 h 1 c 1 i a 2 


—2 

T 

B 2 H 1 D 1 F 1 


y=[ D 2 H i c ii c 2 


. [ D 2 H i d i f i] w 


Figure 9 . State representation for 
concatenated systems . 

The equations defining the two systems are 

• \ 


*1 

= Aixi 

+ 

B 1H1 | 

Zl 

= cixi 

+ 

Diui J 

• 




x„ 

= A„x„ 

+ 

B_u„ 

-2 

2-2 


2-2 j 

^2 

= C 2— 2 

+ 

° 2— 2 , 


System 1 


System 2 


The control vector to the first system is determined by the input to the system and 
the input selection matrix. 

u-| = F-|w 


The equations for the first system can be rewritten as 


xi = A-|X 1 + B-|Fiw 

21 = Cixi + DiF iw 
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The control vector for the second system is defined by the output selection matrix 
and the observation vector from the first system. 


= HiC-jX-) + H i D F -j w 

The state vector x, the control vector u, and the observation vector y for the total 
system are defined as 



Z = Z 2 

The system equations for the second system must be reformulated in terms of xi , 
—2' Z 2 ' an ^ t ' ie system and selection matrices. Using the expanded definition of 
the control vector for the second system, the reformulated equations become 


2 = ^ 2*2 B 2 H iD -jF -jw 

Z 2 = C 2 x 2 + + D 2 H 1 D 1 F •jW 


Using these equations and the expanded equations for the first system (in terms of 
w), the state equations for the total system can be formulated. 



A 1 | 0 1 {"2.1 1 f B 1 F 1 

B 2 HiC-, I a 2 _ _x 2 _ _b 2 h 1 d 1 f 1 _ 


"*r 

_- 2 _ 


Z - [d 2 h 1 c 1 I c 2 ] -- + [d 2 h 1 d 1 f 1 ] w 


Systems With a Common Input 

Figure 10 illustrates two systems with a common input preceded by an input selec- 
tion matrix. The formulation of the partitioned matrices for the total system is 
straightforward once the state equations are rewritten in terms of the input vector 
and input selection matrix. 
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o 


"o 

X 1 


> 

o 

*1 


Vi' 

"o ~ 

x 2 


L°iM 

* 

+ 

B 2 F 1 

V 


N "1 

' 1 

*1 


D 1 F 1~ 

12 


1 

CM 

O 

O 

*2 

T 

D 2 F 1 


Figure 10. State representation 
of systems with a common input. 


The equations for the two systems are. 

X-] = A 1 JC ^ + B iu 
Z1 = C-|Xi + D-iu 


> System 1 


—2 = ^ 2—2 + B 2 — 


Zo = C 2 x 2 + D 0 u 


System 2 


The control vector (u) for the two systems is defined in terms of the input vector w 
and the input selection matrix F-). 


U = F 


1W 


The state vector x, the control vector u, and the observation vector for the total 
system are defined as follows. 


x = 


'*1 


. — 2 . 


u = w 


Z = 


"zr 
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Using the definition of the control vector common to both systems, the system 
equations can be rewritten. 


2£l 

II 

> 

$ 

+ 

B^F^w 

Zl 

= C^x-| 

+ 

D^F ^w 

• 

-2 

= A„x„ 
2-2 

+ 

®2 F 1 — 

—2 

= C.x„ 
2-2 

+ 

D 2 F 1* 


System 1 
System 2 

Thus, by inspection, the total system can be formulated as 


51 


f A 1 1 

i 

°i 


"*r 


_ B 1 F 1~ 

• 

L^J 

— 

I 

1 

0 j a 2 _ 


2_ 

+ 

_ b 2 f 1_ 

~y.il 


o 

O 


5.1 ~ 


_D i F r 

: 

= 

1 

1 


— 

+ 

— 

1 

K 

tsj 

l 


L° 

1 C 2_ 


— 1 
CN 
XI 
1 


) 

c\i 
Q 
I 


Summation of Systems With a Common Input 


Figure 11 illustrates the summation of two systems having a common input. The 
total system is preceded by an input selection matrix. Output selection matrices 
are appended to each system before the summation. After rewriting both sets of 
system equations to include the input vector and input selection matrix, the total 
system formulation follows immediately from the definition of the total system out- 
put vector. 

The systems are defined by the standard equations. 


• 1 

= A-jX-j + BfU 

XI = C-|Xi + DfU 

k 2 = A 2 x 2 + B 2 u 

*2 = C 2 — 2 + ° 2 — , 


System 1 


System 2 
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Figure 11. State representation of summation 
of systems with a common input. 


The common control vector u for two systems is defined in terms of the input vector 
w and the input selection matrix F-j. 

u = F-|w 

The total system output vector can be written in terms of the output vectors of the 
two systems , 

L = £1 + z 2 

or expanded in terms of the output selection matrices and observation vectors of the 
two systems . 

z = H izi + h 2 z 2 

The state vector x and the control vector u for the total system are defined as 
follows . 



u = w 


Rewriting the state equations in terms of the input vector yields 


x-| = A-jx-i + B-|Fiw System 1 
x? = A2X2 + B2F1W System 2 


1 

! 
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The observation vector can be written by substituting for the observation vectors in 
the expanded observation vector definition. 


y — H -)C i jc-| + ^ 2 ^ 2—2 t H 2 ^ > 2 ^' 1 ^ 


Using these equations, the total system can then be formulated in terms of par- 
titioned matrices. 


jy_ 

• 

2_ 

= 

A i | 0 - 

_°”l A 2_ 

"*i; 
. 2_ 



1 

-xi- 


H 1 G ^ 

| «2C 2 ] 

_-2_ 

+ [h 


B 1 F l' 


b 2 f 1 


w 


+ [H^F, + f^F^ 


w 


Feedback System 

Figure 12 illustrates a generalized feedback system with a linear subsystem in 
the forward and feedback loops. An input selection matrix precedes the total system. 
The observation vector of the forward-loop system is used as the total system output 
vector. The feedback system incorporates an input and an output selection matrix. 
After a definition of the total system state, control, and observation vectors, a 
multistep process is used to rewrite the subsystem equations in terms of the elements 
of the total system vectors. 


o 




o 


A^B^jDgFjNC., 

- B 1 H 2 ( C 2- D 2 F 2 ND 1 H 2 C 2 )" 

_ 

£i 



b 2 f 2 nCi 

a 2 - b 2 F 2 nd i h 2 c 2 

- 2 


B 2 F 2 NDiFi 


y = 


NC. 


!-ND 1 H 2 C 2 



w 


Where N = (l+D.,H 2 D 2 F 2 ) 1 

Figure 12. State representation of feedback system. 
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The equations for the 


The control vector to the 
input vector w, the total 

ui = u - 

= F iw - z ^ 2 

This can be expanded in terms of the observation vector of the second system, y_ 2 , 
and the output selection matrix of the second system, H 2 , using the definition 

—2 = V 2 

and substituting into the equation for the control vector for the first system, u-) , 
yields 

ui = f-|w - H 2 y _ 2 

The control vector for the second (feedback) system is defined as 

—2 = Vi 

The state, control, and observation vectors for the total system are defined as 
follows. 


two subsystems are given as 

t ' 

x-) = A-jX 1 + B -|U 

Z1 = + D 1 u 1 J 


/ System 1 


x„ = Ax + B u 
-2 2-2 2-2 


y_ = C x + D u 
X 2 2-2 ">-■ 


2-2 ; 


System 2 


first system, u-| , is defined in terms of the total system 
system selection matrix F-|, and the feedback vector z . 
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The equations for the total system are derived by first expanding the observation 
equation for the forward-loop system using the expanded equation for the control vec- 
tor of the first system. 

XI = c l£l + D l( F l”. “ h 2 x 2 ) 

= Cix-, + D^w - DiH^ 

Substituting the observation equation for the second system into this equation 
yields 


Xl = C 12£1 + D 1 F 1W ~ D-|H 2 (p2— 2 + °2~2) 

= Cix.) + D-|F-|W - D-|H 2 C 2 x 2 - D-|H 2 D 2 u 2 

This equation is expanded again in terms of the definition for the control vector 
for the feedback system 

XI = C-|X1 + D -|F i w - D-|H 2 C 2 x 2 - D-iH^^y.,) 

= C 1 x 1 + DiF-|W - D-|H 2 C 2 x 2 - DiH 2 D 2 F 2 y-, 

which can be written as 

(I + D-|H 2 D 2 F 2 )yi = C-|Xi + D-|F-|W - DiH 2 C 2 X 2 

By defining an intermediate matrix, N, for notational convenience, as 

N = (I + D-|H 2 D 2 F 2 )~ 1 

the observation equation for the forward- loop system can be written 

y 1 = NC ixi - NDiH 2 C 2 X 2 + ND-|F-|W 

In a similar manner, the observation equation for the feedback system can be 
derived. Starting with the simple equation for the observation vector of the second 
system, x 2 > 


X2 = C 2 x 2 + D 2 u 2 


and substituting for u 2 yields 


X 2 = C 2 x 2 + D 2 F 2Z1 
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This equation can be expanded in terms of yq to give 


* 2 = C 2-2 + ° 2 F 2 ( NC 1-1 “ nd 1 h 2 C 2 ^2 + ND 1 F 1 ^) 

= D 2 F 2 NC 1 x 1 + (C 2 - D 2 F 2 ND 1 H 2 C 2 )x 2 + D 2 F 2 ND 1 F 1 w 

Using this equation for the observation vector of the feedback system, the state 
equation for the first system can be written. Using the original state equation 

x i = A-jX^ + B -|U 

and the expanded definition of the control vector for the first system, 

Hi = F iw - H 2 y _ 2 

the state equation becomes 

xi = A-|x-| + B -| ( F iw - H 2 y 2 ) 

= A-jxi + B-|F-|W - BiH 2 y_ 2 

which can be expanded in terms of the observation vector for the second 
system, ^ , to give 

xi = Aix-| + B-|F -|w 

- B i H 2 [D2F 2 NC-|Xi + (C 2 - D 2 F2ND-,H2C2 )x 2 + D 2 F 2 ND 1 F 1 wJ 
Rearranging terms yields 

xi = ( A-| - BiH 2 D 2 F 2 NC'j )xi - B-|H 2 (C 2 - D 2 F 2 ND -|H 2 C 2 )x 2 
+ (B -|F 1 - BiH 2 D 2 F 2 ND-|Fi )w 

Finally, an expression for the state equation for the second system is derived. 
This derivation is based on the definition of the control vector for the feedback 
system u 2 , 

H 2 = F 2 Z 1 

and the final expression for yq , 

yi = NC 1X1 - NDiH 2 C 2 x 2 + ND -)F -|W 
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which expands the definition to 


\i — F 2NC ^ x*| “ F2ND 1H2C x + F 2 ND-jF<jW 
^ 4 z 

This expression for u^ is substituted into the original state equation for the feed- 
back system 


—2 


= Ax 
2-2 


+ B u 
2-2 


to yield 


* 2 = A 2 -2 + B 2 ( F 2 NC 1*1 " f 2 nd 1 h 2 c 2 * 2 + F 2 ND 1 F 1 ”) 

Rearranging terms results in a usable expression for the state equation of the 
second system. 


x 2 = B 2 F 2 NC 1 x 1 + 


( a 2 ~ B 2 F 2 ND 1 H 2 c 2)2i 2 + B 2 F 2 ND 1 FiW 


The final state equations and the observation equation for the first system can 
be expressed in terms of partitioned matrices for the total system equation. 


• “ 

-1 


fA-i - B-|H 2 D 2 F 2 NCi 

1 -BiH 2 (c 2 - d 2 f 2 nd -|H 2 C 2 ) ”1 


"*r 

— — 




— 



1 b 2 f 2 nci 

1 A 2 _ B 2 F 2 ND 1 H 2 C 2 J 


L— 2 J 


+ 


B 1 F 1 - BiH2D2F 2 NDiFi 
b 2F 2 NDiFi 


w 


Z = [ NC 1 J -NDiH^] 



+ f ND -)F -j] w 


MATRIX x COMMAND FILE 

To aid in the process of combining and connecting linear multi-input, multi- 
output subsystems into a total system model, a command file has been created in 
MATRIX x format. This file uses subsystem matrices whose names are consistent with 
the notation of this report and produces total systems matrices ANEW, BNEW, CNEW, 
and DNEW which correspond to the A, B, C, and D matrices, respectively. The options 
provided by this command file are the six sample systems developed in this report. 

In addition to the four matrices above, the command file produces a total system 
matrix, SNEW, defined in the following way. 
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SKEW - 


“ANEW I SNEW‘ 
JcNEW* "t" "dNEWL 


While the command file combines subsystems two at a time, the file is designed to 
allow multiple combinations without leaving the command file. 

The command file is used within MATRIX x by typing 

EXEC ( ' COMBLK ' ) 


The command file responds by typing a title, 

COMMAND FILE TO ASSIST IN COMBINING AND CONNECTING SYSTEMS 
REFERENCE NASA TM-85912 


After the title is printed, the screen is cleared and a menu is displayed. 

WANT TO DEFINE MATRICES 

( 1 ) NO 

(2) YES 

The selection is made by positioning the cursor next to the desired answer and then 
depressing the RETURN key. 

If all matrices have been named using the notation of this report, the 
matrices do not have to be redefined. However, if other names have been 
used, the command file will request that the user define the appropriate 
matrices. A typical response following an answer of "yes" to the question 
"want to define blocks" is 

ENTER THE STATE MATRIX FOR THE FIRST SYSTEM, A1 

The user must then either enter the elements of the A^ matrix or equate the Ai 
matrix to another, previously entered matrix, representing the state matrix for the 
first system. 

After the matrices for both systems are defined in the correct notation, the 
screen is again cleared, and another menu appears. 

SELECT AN OPTION 

( 1 ) PARALLEL SYSTEMS 

(2) SUMMATION 

( 3 ) CONCATENATION 

(4) COMMON INPUT SYSTEMS 

(5) SUM COMMON INPUT SYST 

( 6 ) FEEDBACK SYSTEM 
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Again, the selection is made by positioning the cursor next to the desired option 
and then depressing the RETURN key. The options listed correspond to the examples 
detailed in this report. 

If the user has elected to enter or define the matrices, the command file will 
again request that the required input and output selection matrices be entered. The 
actual matrices to be defined are determined by the option selected. However, 
because all options require an input selection matrix for the first system, F-j, the 

following message will always appear if the user is entering or defining matrices. 
ENTER THE INPUT SELECTION MATRIX FOR THE FIRST SYSTEM, FI 

The user must then either enter the elements of the F-| matrix or equate the Fi 

matrix to another, already entered matrix, representing the input selection matrix 
for the first system. 

The screen will again be cleared and another menu will appear. 

DO YOU WANT TO EXIT 

( 1 ) NO 
( 2 ) YES 

The selection is made by positioning the cursor next to the desired answer and 
depressing the RETURN key. If (1) is selected, the command file will return to the 
first menu and again ask if the user desires to define the matrices. Selection of 
(2) will cause the system to exit the command file and return the user to MATRIX x . 


Ames Research Center 

Dryden Flight Research Facility 

National Aeronautics and Space Administration 

Edwards, Calif., August 2, 1984 
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APPENDIX A - LISTING OF COMMAND FILES 


There are eleven command files to assist in combining and connecting systems. 
These command files each perform specific functions described below. The prolifera- 
tion of these files is partly due to normal program structuring and partly the 
result of a 1024-character limitation on command files within MATRIX x . 

The COMBLK command file (fig. A. 1 ) is the main file used to combine and connect 
systems. This file acts as an executive and controls the basic operation of the 
command file. MENUS (fig. A. 2) is a command file to define the menus listed when 
using the COMBLK command file. The INPSYS command file (fig. A. 3) is used to pro- 
vide the user prompting while entering the system matrices. The input of selection 
matrices is controlled by the SELECT command file (fig. A. 4). SELECT contains logic 
to request only the selection matrices required for the application selected by the 
user. The actual combining and connecting of systems is controlled by the EXCOPT 
command file (fig. A. 5). 

Figures A. 6 to A. 11 list the command files used to determine the matrices for 
the total combined system. The combining of parallel . systems into a single system 
description is done by the PARSYS command file (fig. A. 6). The summation of 
parallel systems is performed by the SUMSYS command file (fig. A. 7). The CONSYS 
command file (fig. A. 8) determines the total system matrices when the systems are 
to be concatenated. The COMINP (fig. A. 9) and SUMCOM (fig. A. 10) command files 
determine the total system matrices for systems with a common input and separate 
outputs or with a common input and a summed output. The FEDBAK command file 
(fig. A. 11) determines the total system matrices for a feedback system. 
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// 


// COMMAND FILE TO ASSIST IN COMBINING SYSTEMS 

// 

// REFERENCE NASA TM-85912 

// LEE DUKE 

// 

EXEC ( 'MENUS' ) 

DISP( 'COMMAND FILE TO ASSIST IN COMBINING AND CONNECTING SYSTEMS') 
DISP( 'REFERENCE NASA TM-859 1 2 ' ) 

ANSWER3=1 

WHILE ANSWR3=1 ; . . . 

ANSWR1=MENU (LISTOP ) - 1 ; . . . 

IF ANSWR1=1 ;EXEC( ' INPSYS ' ) ; END ; . . . 

ANSWR2=MENU ( LISTCM ) ; . . . 

IF ANSWR1=1;EXEC( 'SELECT' ) ;END; . . . 

EXEC( 'EXCOPT' ) ; . . . 

ANSWR3=MENU( LISTCH ) ; . . . 

END; 

Figure A.l Listing of COMBLK command file. 
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// 


// 

// 

LISTOP 


LISTCM 


LISTCH 


COMMAND FILE TO SET UP MENUS FOR COMBINING SYSTEMS 

= ['WANT TO DEFINE BLOCKS' 

'NO ' 

'YES '] 

= ['SELECT AN OPTION ' 

'PARALLEL SYSTEMS ' 

' SUMMATION ' 

' CONCATENATION ' 

'COMMON INPUT SYSTEMS ' 

'SUM COMMON INPUT SYST ' 

'FEEDBACK SYSTEM ']; 

= ['DO YOU WANT TO EXIT' 

'NO ' 

'YES ']; 


Figure A. 2 Listing of MENUS command file 


// 


// COMMAND FILE TO INPUT SYSTEM MATRICES 

// 

DISP ('ENTER THE STATE MATRIX FOR THE FIRST SYSTEM, A1 ' ) ; 

EXEC ( 5 ) ; 

DISP ('ENTER THE CONTROL MATRIX FOR THE FIRST SYSTEM, B 1 ' ) ? 

EXEC ( 5 ) ; 

DISP ('ENTER THE OBSERVATION MATRIX FOR THE FIRST SYSTEM, Cl') 
EXEC ( 5 ) ; 

DISP ('ENTER THE FEEDFORWARD MATRIX FOR THE FIRST SYSTEM, Dl')? 
EXEC ( 5 ) ; 

DISP ('ENTER THE STATE MATRIX FOR THE SECOND SYSTEM, A2 ' ) ; 

EXEC ( 5 ) ; 

DISP ('ENTER THE CONTROL MATRIX FOR THE SECOND SYSTEM, B2 ' ) ; 
EXEC( 5 ) ; 

DISP ('ENTER THE OBSERVATION MATRIX FOR THE SECOND SYSTEM, C2 ' ) ; 
EXEC ( 5 ) ; 

DISP ('ENTER THE FEEDFORWARD MATRIX FOR THE SECOND SYSTEM, D2 ' ) ; 
EXEC ( 5 ) ; 

Figure A. 3 Listing of INPSYS command file . 
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// 


// COMMAND FILE TO INPUT SELECTION MATRICES 

// 

DISP ('ENTER THE INPUT SELECTION MATRIX FOR THE FIRST'); 

DISP (’SYSTEM, FI'); 

EXEC ( 5 ) ; 

IF ANSWR2 < > 3; IF ANSWR2 < > 4; IF ANSWR2 < > 5; . . . 

DISP ('ENTER THE INPUT SELECTION MATRIX FOR THE SECOND'); . . . 

DISP ('SYSTEM, F2 ' ) ; . . . 

EXEC ( 5 ) ; . . . 

END; END; END; 

IF ANSWR2 < > 1; IF ANSWR2 < > 4; IF ANSWR2 < > 5; . . . 

DISP ('ENTER THE SELECTION MATRIX FOR THE OUTPUT OF THE FIRST');. . 
DISP ('SYSTEM, HI'); . . . 

EXEC ( 5 ) ; . . . 

END; END; END; 

IF ANSWR2 < > 1; IF ANSWR2 >4; . . . 

DISP ('ENTER THE SELECTION MATRIX FOR THE OUTPUT OF THE SECOND');. 
DISP ('SYSTEM, H2 ' ) ; . . . 

EXEC( 5 ) ; . . . 

END ; END; 


Figure A. 4 Listing of SELECT command file 



// 

// COMMAND FILE TO EXECUTE OPTIONS FOR COMBINING 
// AND CONNECTING SYSTEMS 
// 


IF 

ANSWR2 = 1? 

EXEC 

( 'PARSYS' 

) ; 

END? 

IF 

ANSWR2 = 2; 

EXEC 

{ 'SUMSYS' 

) ; 

END; 

IF 

ANSWR2 = 3; 

EXEC 

( 'CONSYS ' 

) ; 

END? 

IF 

ANSWR2 = 4; 

EXEC 

( ' COMINP ' 

) ; 

END; 

IF 

ANSWR2 = 5; 

EXEC 

( ' SUMCOM ' 

) ; 

END; 

IF 

ANSWR2 = 6; 

EXEC 

( 'FEDBAK' 

); 

END? 


SNEW = [ANEW BNEW; CNEW DNEW] ; 


Figure A. 5 Listing of EXCOPT command file. 


// 

/ / COMMAND FILE TO COMBINE PARALLEL SYSTEMS 

// 

<N1 ,M1 > = SIZE(BI); <Rl,Ml> = SIZE(DI); <M1,L1> = SIZE(FI); 
<N2,M2> = SIZE ( B2 ) ; <R2,M2> = SIZE(D2); <M2,L2> = SIZE(F2); 
ANEW = [A1 0*ONES(N1,N2) ; 0*ONES ( N2 , N 1 ) A2] ? 

BNEW = [B 1 *F 1 0*ONES(N1 ,L2) ; 0*ONES ( N2 , L 1 ) B2*F2] ; 

CNEW = [Cl 0*ONES(R1 ,N2) ; 0*ONES ( R2 ,N 1 ) C2] ; 

DNEW = [D 1 *F 1 0*ONES ( R 1 , L2 ) ; 0*ONES(Rl ,Ll ) D2*F2]; 


Figure A. 6 Listing of PARSYS command file. 
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// 

// COMMAND FILE TO SUM PARALLEL SYSTEMS 

// 

<N 1 / Ml > = SIZE ( B 1 ) ; <M1,L1> = SIZE(FI); 

<N 2 , M2 > = SIZE(B2); <M2,L2> = SIZE(F2); 

ANEW = [A 1 0*ONES(N1,N2) ; 0*ONES (N2 ,N 1 ) A2] ; 

BNEW = [B 1 *F 1 0*ONES(N2,L2) ; 0*ONES (N2 , Li ) B2*F2]; 
CNEW = [H 1 *C 1 H2*C2 ] ; 

DNEW = [H 1 *D 1 *F 1 H2*D2*F2 ] ; 

Figure A. 7 Listing of SUMSYS command file. 

// 

// COMMAND FILE TO CONCATENATE SYSTEMS 

// 

<N 1 / M 1 > = SIZE ( B 1 ) ; 

<N2,M2> = SIZE(B2 ) ; 

ANEW = [A2 0*ONES(N1,N2) ; B2*H1*C1 A2] ; 

BNEW = [B 1 *F 1 ; B2*H 1 *D 1 *F 1 ] ; 

CNEW = [D2*H 1 *C 1 C2] ; 

DNEW = [D2*H 1 *D 1 *F 1 ] ; 

Figure A. 8 Listing of CON SYS command file. 
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// 

// COMMAND FILE TO COMBINE FEEDBACK SYSTEM 

// 

<N1 ,Ml> = SIZE(Bl) ; 

N = INV ( EYE( N 1 ) + Dl*H2*D2*F2) ; 

ANEW = FaI-B 1 *H2*D2*F2*N*C 1 -B 1 *H2 * (C2-D2*F2*N*D 1 *H2*C2 ) • 
B2*F2*N*C 1 A2-B2*F2*N*D1*H2*C2] ; 

BNEW = [B1*F1-B2*H2*D2*F2*N*D1*F1 ; B2*F2*N*D 1 *F 1 ] ; 

CNEW = [N*C1 -N*D1*H2*C2] ; 

DNEW = [N*D 1 *F 1 ] ; 


Figure A. 11 Listing of FEDBAK command file 
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